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ABSTRACT 



Aims. We test the isotropy of the Hubble diagram. At small redshifts, this is possible without assumptions on the cosmic inventory and provides 
a fundamental test of the cosmological principle. At higher redshift we check for the self-consistency of the ACDM model. 
Methods. At small redshifts, we use public supemovae (SNe) la data to determine the deceleration parameter qo and the SN calibration on 
opposite hemispheres. For the complete data sets we fit SIm and the SN calibration on opposite hemispheres. 

Results. A statistically significant anisotropy of the Hubble diagram at redshifts z < 0.2 is discovered (> 95%C.L.). While data from the 
North Galactic hemisphere favour the accelerated expansion of the Universe, data from the South Galactic hemisphere are not conclusive. The 
hemispheric asymmetry is maximal toward a direction close to the equatorial poles. The discrepancy between the equatorial North and South 
hemispheres shows up in the SN calibration. For the ACDM model fitted to all available SNe, we find the same asymmetry. 
Conclusions. The alignment of discrepancies between hemispheric Hubble diagrams with the equatorial frame seems to point toward a sys- 
tematic error in the SN search, observation, analysis or data reduction. We also find that our model independent test cannot exclude the case of 
the deceleration of the expansion at a statistically significant level. 

Key words, observational cosmology - large scale structure - supemovae 
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1. Introduction 

The Hubble law is a direct consequence of the cosmological 
principle. Modern Hubble diagrams from supemovae la (SNe 
la) confirm the Hubble law and provide evidence for an accel- 
erated expansion of the Universe (Riess et al.[T998; Perlmutter 
et al. |1999t . In these studies the isotropy of the Hubble diagram 
is assumed. The purpose of this work is to provide quantitative 
tests of the isotropy of SNe la Hubble diagrams, beyond the 
identification of the cosmic microwave dipole in the local SNe 
la data (Riess et al. U995b . 

A test of the isotropy of Hubble diagrams is interesting 
' for at least three reasons: 1. for checking the validity of the 
cosmological principle (Kolatt & Lahav 2001 ), 2. to measure 
expected deviations from the isotropy [due to Local Group 
motion, peculiar velocities, and other effects from structure 
formation (Sasaki IT9871 Radburn-Smith et al. [20041 Bonvin 
at al. I2006al Hui & Greene [20061 Weinhorst [20061 Cooray 
& Caldwell [20061 Haugb0lle et al. [20061 Neill et al. 120071 
Wang 120071 Hannestad et al. 120071 Gordon et al. 120071) 1. 3. to 
search for systematic errors in the observations and their anal- 
ysis (Kolatt & Lahav I200T1 Gupta et al . 120071 . 

The cosmological principle states that the statistical distri- 
bution of matter and light in the Universe is isotropic and ho- 
mogeneous in space. The principle itself may either be moti- 
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vated by the idea of cosmological inflation or by a simplicity 
argument. As a consequence of the ergodic theorem (spatial av- 
eraging replaces ensemble averaging) the spatial isotropy and 
homogeneity of the statistical distributions becomes an approx- 
imate symmetry of the space-time metric and matter distribu- 
tion at large scales. 

Consequently, the cosmological principle implies that dis- 
tances, angles, time intervals etc. at large scales can be mea- 
sured according to the ruler sticks and clocks described by the 
Robertson- Walker line element. The redshift of light from dis- 
tant (with the expansion comoving) objects and the Hubble 
law for redshift z <k 1 follow directly (without making use 
of Einstein's equation), i.e. the Hubble law does not depend on 
the details of the cosmological model (like the inventory of the 
Universe). The Hubble diagram (distance or magnitude versus 
redshift) therefore provides one of the most fundamental tests 
of modern cosmology. It provides a test of the cosmological 
principle and the idea that the space-time is correctly modelled 
as a riemannian manifold. 

However, the local Universe is neither isotropic nor homo- 
geneous. Consequently, the observations are expected to ap- 
proach the Hubble law at some z > only. Deviations from 
the isotropic and homogeneous Universe arise from structure 
formation and can be described by linear perturbation theory 
at large scales. Thus in order to measure the Hubble constant, 
the deceleration parameter and other cosmological parameters 
of an isotropic and homogenous model, it is compulsory to 
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demonstrate that the effects from structure formation are either 
negligible or properly considered in the error bars. 

Most importantly, one would like to firmly establish the ac- 
celeration of the Universe in a model independent way; i.e. to 
measure the kinematics of the Universe without assuming the 
validity of Einstein's equation and without assumptions on the 
matter content of the Universe. This leads to a dilemma: at very 
small redshifts, the Hubble law does not hold, since the local 
Universe is inhomogeneous. At high redshift, z ~ 1, the lumi- 
nosity distance c/l (or equivalently the distance modulus m-M) 
depends on the detailed cosmological model. Therefore, such 
an analysis must be restricted to a finite range of redshifts be- 
low one, but must exclude the most local SNe. 

A simple and powerful method to establish the existence of 
such a range of redshifts is to test the isotropy of the Hubble di- 
agram. To be more precise, the isotropy of the Hubble diagram 
is a necessary, but not a sufficient condition for the existence of 
such an interval. Here we use SNe la data to test the isotropy 
of the Hubble diagram on opposite hemispheres of the sky. 

This test is also closely related to the issue of fairly sam- 
pling the Universe with SNe. In other words, is it suitable to 
determine the global Hubble rate from a local measurement. 
The 1998 cosmology-revolution (Riess et al. 119981 Perlmutter 
et al. 1999 ) relies on the assumption that the SNe la at low red- 
shifts (say z < 0.2) represent a fair sample of the Universe and 
that the Hubble law is a good approximation to the data. Based 
on SNe la at low and high redshifts together, it has been con- 
cluded that the present expansion of the Universe accelerates. 

At the same time, our test provides a cross check for the 
measurements of the Hubble constant Ha and the deceleration 
parameter qo- The most commonly adopted values for the infla- 
tionary A cold dark matter model (ACDM) are from the HST 
key project Hq = 72 + 8 km/s/Mpc (Freedman et al. 2001) 
and from a fit to WMAP data H = 73^ km/s/Mpc and 
q = (3D.M - 2)/2 = -0.64+° °^ (derived, not directly mea- 
sured) (Spergel et al. (2006). The calibration of SNe la used 
in the HST key project analysis has recently been criticised by 
Sandage et al. 2006 The HST key project analysis is based on a 
set of cepheides from a metal-poor environment, whereas most 
of the SNe of the sample are in metal-rich galaxies. Sandage et 
al. 12006 obtain Hq = 62. 3± 1.5 ±5.0 km/s/Mpc, using cepheides 
from a metal-rich environment and a significantly larger num- 
ber of calibrators and SNe. At face value, there is now some 
tension between CMB measurements and the direct determina- 
tion from the Hubble diagram. A possible explanation would 
be that the local value does not coincide with the global one. 
Our test is suited to detect such a situation, as we would expect 
it to go along with an anisotropy of the Hubble diagram. 

In the second part of this work we present a model depen- 
dent analysis of the isotropy of Hubble diagrams, relying on 
the spatially flat ACDM model, i.e. a model with two free pa- 
rameters Ho and Qm- 

The paper is structured as follows. In section|2]we comment 
on some basic properties of the Hubble law. In section [3] we 
describe the four SNe la data sets that we are using. The method 
of our tests is explained in section |4] and section [5] contains the 
results. We conclude in section [6] 



2. Hubble law 

The luminosity distance d^ is a function of redshift and angu- 
lar position of an observed object. It is linked to the distance 
modulus m - M — 5 log^/ IMpc) + 25. The observed redshift 
z of the object contains information on its peculiar velocity, the 
peculiar velocity of the observer and the cosmological redshift. 
If we assume that the cosmological principle holds and that 
the observed objects are comoving with the expansion of the 
Universe, one finds the Hubble law 
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where j denotes the jerk, k = 0, + 1 the normalised, comoving 
curvature and a the scale factor. In terms of the scale factor and 
its derivatives 
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We see that the luminosity distance does not depend on the 
geometry of the Universe up to order z 2 , while the terms of 
order z 3 depend on jerk and geometry. We restrict our analysis 
to a test of the Hubble law up to order z 1 ■ A useful reference 
is the empty-Universe model, which is characterised by q — 
j = 0. Using the Einstein equation for the reference model we 
find additionally H 2 /c 2 = -k/a 2 , thus k = -1 and Hod^iz) = 
cz(l+z/2). 

Let us estimate the error from neglecting jerk and curvature 
in ([T). For a general ACDM model, we find qo = (Qm -2Qa)/2 
and jo = £2m + Qa- In the special case of the Einstein-de Sitter 
model (qo = 1/2, jo = l,k - 0), the ratio of the third order 
term to the second order term in ([]]) becomes -1/2 and thus 
the error from neglecting the third order term amounts to 10% 
at z = 0.2. For the WMAP best-fit flat ACDM model the error 
at the same redshift is 6%. Therefore, in order to obtain qo and 
its confidence contours at 10% theoretical accuracy or better, 
we must restrict the model independent fits to SNe at redshifts 
below z = 0.2. 

For the spatially flat ACDM model (assuming Einstein's 
equation) the exact expression for the luminosity distance is 
given by 



d h (z) = c(l + z) 



dz' 
Hiz'Y 



(3) 



with H(zf = H 2 [D. W (\ + z) 3 + 1 - Q M ]- Thus this model is 
fully characterised by the parameters Ho and Qm, which can be 
easily compared to the model independent parameters Hq and 
qo- Expression ^ is the basis of our model dependent fit and, 
in contrast to ([TJ, holds for arbitrary redshift. 

The effect of peculiar velocities on the redshift can be in- 
corporated by realising that 

c(l + Zobs) = c(l + z com ) [l + e(v pec - v )/c + ...], (4) 

where e denotes the direction of observation (unit vector) — 
light propagates along -e — and the dots stand for higher or- 
der terms in the peculiar velocities. As shown by Bonvin et 
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al. d2006al I2006bl l. it is expected that the peculiar motion of 
the Solar system w.r.t. the CMB is the dominant source of 
anisotropy in the Hubble diagram. Indeed the CMB dipole has 
its correspondence in SN la data (Riess et al. 119951 1. It is eas- 
ily subtracted, assuming that the CMB dipole is entirely due 
to the Sun's motion w.r.t. the CMB. We take the WMAP lyr 
value for the dipole: (l,b) = (263.85 o +0.r,48.25 o +0.04°) and 
v Q /c = Ti/Tq = [(3. 346+0.017)/(2.725±0.001)] = 1.23xl(T 3 
(Bennett et al. 2003 ). At small redshifts it is sufficient to work 
with 

CZcmb ~ CZcom t-Vp e c + • ■ • (5) 

We assume that the peculiar velocities are uncorrected with 
zero mean and variance crj. 

For the purpose of our model independent analysis we de- 
fine the deviation from the linear Hubble law 

A = \og(H Q d L ) - log(cz cmb ), (6) 

where log = log 10 . In the following z cm b is denoted by z. Note 
that often the empty Universe model is plotted as a reference, 
but that model assumes already the Einstein equation and we 
would like to avoid in our first test any reference to the dynam- 
ics of the Universe and stick to a purely kinematic test com- 
bined with the spatial information provided by the data. 

The quantity A is free of the Hubble constant, once (Q} or 
(0 are used. Equivalently, log(//o^L) is related to the "Hubble- 
constant-free" distance modulus 

m - M = 5 log(H d L ), 

with M = M - 5 log H + 25 (Perlmutter et al. 1 19991 . However, 
to specify the SN data (in terms of log(//o<i) or m - M) a cer- 
tain calibration has to be adopted, which makes use of fiducial 
values Hq and M*. 

In order to use SNe for the study of cosmology, one as- 
sumes that type la SNe can be made standard candles to a good 
approximation, i.e. Al becomes a universal number. As this 
number depends mainly on the physics of SNe (the universality 
of Hq is taken for being granted), it is treated as a nuissance pa- 
rameter in the cosmological data analysis. Consequently, for 
the estimation of cosmological parameters one marginalises 
over M (or equivalently Hq). However, if we are interested in 
effects of large-scale structure formation or a test of the cos- 
mological principle, we need to study the direction dependent 
off-set from the adopted calibration. 

Below we follow the convention of Tonry et al. (2003), 
where a Hubble rate of Hq = 65 km/s/Mpc is assumed; we 
thus define (the index i denotes a particular SN) 

Ai = log(#^ L i)-log(cZi), (7) 

which is compared to the second order expression 
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Fig. 1. Distribution of SNe la on the sky in galactic coordi- 
nates. Red squares denote SNe with z < 0.2, whereas blue disks 
stand for SNe with z > 0.2. Only SNe with extinction A v < 1 
are shown. For SNe indicated by open symbols no information 
on Ay is available. From top to bottom: 244 (105 with z < 0.2, 
139 with z > 0.2) SNe la from Tonry et al. (l20"03l ) and Barris et 
al. (120041) : 182 S Ne la ( 77 with z < 0.2 and 105 with z > 0.2) 
from Riess et al. (l2QQ4t ; 1 17 SNe la (44 with z < 0.2 and 73 
with z > 0.2) from Astier et al. (2006), one can easily spot 
the four SNLS fields of view (blue); 119 SNe la from Jha et 
al. (120061 ) all at small redshift. 



A*( Zi ) = - log(H /H* ) + l °SU + (1 - «o)f L (8) H o/H*, or M - M* = -5 \og(H /H*), i.e. a calibration off- 

set. Thus, a 10% effect in Ho /HI corresponds to 0.2 mag in M. 
or the equivalent in case of the flat ACDM fits. We treat all data The difference M - AT seems to be better suited for the corn- 
sets (see below) in the same way. parison of different data sets, whereas we work with Hy/H*, 
As we have no access to the absolute calibration of the SNe, which seems to us better suited for the study of cosmological 
we cannot measure Hq. However, we can measure the ratio anisotropies. 



4 



Dominik J. Schwarz and Bastian Weinhorst: (An)isotropy of the Hubble diagram: comparing hemispheres 




set reference(s) 



full set 
#SNe (z) 



A V <1 
#SNe (2) 



Tonry et al. 03 
Barris et al. 04 

B Riessetal. 04 186(0.393) 182 (0.398) 77 (0.041) 

C Astier et al. 06 117(0.407) 117(0.407) 44(0.038) 

D Jhaetal. 06 131 (0.023) 119 (0.024) 119 (0.024) 

Table 1. Compilation of some characteristics of data sets and 
their subsamples used throughout this work. 



A v <l,z<0.2 
#SNe (2) 



253 (0.263) 244 (0.273) 139 (0.037) 




Fig. 2. Deviation from the linear Hubble law (see equation |7]i 
for the four data sets from figureQ] SNe from the North (South) 
galactic hemisphere are shown in blue (red). As a reference we 
show the expected deviation for the empty Universe model (q = 
j = Q;k = -1) with three different values of the SN calibra- 
tion (black solid lines), corresponding to Hq/Hq = 0.85, 1, 1.15 
from top to bottom. 

3. Data 

We apply our tests on four different data sets that differ in sky 
and redshift coverage, shown in figures Q] and [2] as well as in 
systematics. A summary of the four data sets is given in table 
ffl 

Our first data set (A) consists of 253 SNe with redshifts 
z e [0.0023, 1.755] (mean z = 0.27) from Tonry et al. (120031 
and Barris et al. ( 120041 1. This data set has a fairly homogeneous 
sky coverage (except in the zone of avoidance) and is therefore 
suitable for our analysis. For each SN la we use the (galactic) 



coordinates, the distance modulus and the redshift in the CMB 
frame. We also use information on light extinction to exclude 
supernovae with Ay > 1.0 from our analysis. A subset of 42 
SNe (originally from the Supernova Cosmology Project) has 
no information on Ay specified. We include the specified errors 
of the distance modulus and the redshift in our analysis and 
treat them as lcr errors. The errors on the distance modulus in- 
cludes the error from photometry, the scatter between different 
analysis methods and an error from the intrinsic dispersion of 
SNe la. 

The second data set (B) is from Riess et al. (120041 1. con- 
taining 186 SNe with z e [0.0104, 1.755] with mean z = 0.40. 
Data set B is not independent from data set A, as it largely con- 
tains the same SNe. However, different corrections have been 
applied and different selection criteria have been used. A com- 
parison of data set A and data set B can thus serve as an es- 
timate of systematic errors in the SNe la data reduction. Data 
set B contains a 'gold' and a 'silver' set. We make use of the 
'silver' set here in order to have a good sky coverage. It should 
be stressed that the number of SNe with z < 0.2 in set B is 
significantly smaller than in set A. 

As a third data set (C) we use the more recent SuperNova 
Legacy Survey (SNLS) release (Astier et al. l2"0"06i l; 117 SNe, 
z = [0.015, 1.01] and z = 0.41. The SNe at high redshift are 
independent from sets A and B, but come from four fields of 
view only (see figure [TJ. Thus we should expect that the pen- 
cil beam geometry of this data set is not ideal for the type of 
analysis that we have in mind, although it provides improved 
control over observational systematics and has very small er- 
rors compared to data set A and B. At low redshifts, data set C 
relies largely on the same SNe as data set A and B, but again 
with different processing and selection criteria, so it provides 
again a cross check for systematics. In contrast to the sets A 
and B, the intrinsic scatter of SNe is not included in the error 
provided by the SNLS team, rather it has to be incorporated as 
an additional term in the fits. 

Finally, we also use the most recent compilation of 'lo- 
cal' SNe from Jha et al. d2006t (data set D); 131 SNe, z e 
[0.00226,0.1242] and z = 0.024. It largely contains the same 
SNe as data set A, the main difference being a largely improved 
procedure to obtain the distance moduli and their errors. This 
set has a fairly good sky coverage, but is limited to small red- 
shifts. Like for data set C the intrinsic scatter is not included in 
the provided error. 

Two more recent data sets are not used because they contain 
SNe from very small regions of the sky only, i.e. from the two 
HST GOODS fields (Riess et al. i2006T > and 20 fields of view 
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Fig. 3. Confidence contours for a model-independent full-sky 
fit to the Hubble law at second order for three SNe la data sets. 
SNe up to redshift z - 0.2 are included in the fits. We determine 
the calibration off-set Hq/H^ and the deceleration parameter 
qo. As a comparison we show the corresponding contours of 
the WMAP measurement (Spergel et al. I2006I I. In contrast to 
the spread of the SN calibration off-set, its best-fit value has no 
physical relevance in the case of full sky fits. Thus the WMAP 
data may be moved horizontally in this figure. 



within three hours in right ascension (ESSENCE collaboration: 
Wood-Vasey et al. 2007). These surveys suffer from the same 
limitations as the SNLS data, they are a set of pencil beams. 

For all data sets we can see from figure Q] that the sky cov- 
erage is much better for the low redshift SN (z < 0.2). Figure 2 
provides an impression of the redshift distribution of the sam- 
ples for both galactic hemispheres. The current lack of SN ob- 
servations at intermediate redshifts (z = 0.1 to 0.3) is easily 
spotted. At small z < 0.01 the spreading of the distribution due 
to peculiar velocities of the SNe is seen. The deviation from the 
linear Hubble law becomes visible at z > 0.1. 

The four data sets used here differ from each other in the 
method applied to estimate the values of the fitted parameters. 
In order to fit the SN observations of set A, Tonry et al. {2003) 
and Barris et al. (120041 used four fit methods: the MLCS 
method described in Riess et al. d 1998) , the Amis method 
(Philips et al. 1999) and its improved version dml5 (Germany 
et al. 1200 11 1. Additionally, they used the BATM method by 
Tonry et al. (2003). For data set B, the BATM method was used 
to check the results with the improved version MLCS2k2 de- 
scribed in Jha (2002) and Jha et al. d20061 l. This method is also 
used by Jha et al. (120061) to fit the SNe of data set D. Astier et 
al. (|2006 ) used the SALT method (Guy et al. |2"0"05T l to fit the 
data of set C. 

The reader should be aware that Tonry et al. (2003 ), Barris 
et al. (120041 ) and Riess et al. ( 120041 ) provide z C mb for each SN, 
whereas Astier et al. (2006 ) and Jha et al. (|2006) give z . 



4. Comparing hemispheres 

As motivated in the introduction, we wish to test the isotropy 
of the Hubble diagram. As the number of observed SNe is still 
small, it does not make sense to look into small-scale varia- 
tions, rather we should look at the largest possible scales. The 
largest possible anisotropy scale is 180°. A very simple test 
is splitting the sky into hemispheres and to compare the cor- 
responding Hubble diagrams. We can search for directions of 
anisotropies by rotating the respective poles over the sky. A 
similar study has been used by Eriksen et al. (2004) to anal- 
yse the cosmic microwave sky, which led to one of the several 
anomalies discovered in the WMAP data (see Copi et al. 2007 
for a recent update and summary). 

4.1. Model-independent test of isotropy 

If we restrict the analysis to SNe at z < 0.2, it is reasonable to fit 
the Hubble law (HJ and its quadratic correction. This provides 
a model independent test of the cosmological principle and at 
the same time a check for systematic errors in the SN la search, 
observation and analysis. 

In order to find the best-fitting values of the calibration 
Ho/H* Q and the deceleration go, we follow the standard ap- 
proach and minimise 



X 2 (Ho/H* , qo ) = J] 



[A / -A*fc;// ///o'4o)f 



2 ' 



(9) 



where the sum goes over a subset of one of the three data sets. 
For the data sets A and B the intrinsic dispersion ctm is in- 
cluded in <TA t , for data set C and D cr m i > 0. cr v and cr, denote 
the dispersion of the peculiar velocities and error of the redshift 
measurement. The latter two are most relevant for nearby SNe. 

4.1.1. Full Sky 

As a reference we first determine the best fit values for Hq/HZ 
and qo for the four SN la data sets. In doing so we aim at finding 
the best suited values for redshift and light extiction cuts, and 
we look a the best suited value of the peculiar velocity disper- 
sion. We adopt a value of cr v = 345 km/s for all four data sets, 
as it provides a reasonable x 2 P er degree of freedom of order 
one for all data sets under consideration. The z < 0.2 samples 
from sets A to D have mean redshifts of 0.038, 0.041 , 0.038 and 
0.024, respectively. 

The results of the full-sky fits are presented in figure[3](ex- 
cept for set C) and in table [2] The fits shown in the figure cor- 
respond to the first fit in the table for the respective data set. 
We find that all four data sets provide fits that are consistent 
with each other. The largest difference is observed with respect 
to the SN calibration (Ho/H^) in data set C. This happens be- 
cause the SN calibration of that set has been fit to the WMAP 
measurement of Ho, whereas the other three data sets did not 
assume the WMAP value. At small redshifts, only data set B 
gives rise to a marginal evidence for acceleration. However, 
all data sets are consistent with acceleration and with the fit to 
WMAP data. 
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data set A data set 
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Fig. 4. North (full lines) and South (dashed lines) confidence contours and best-fit values for galactic, equatorial and maximum 
asymmetry hemispheres for the model-independent test. These fits should be compared to the full-sky fits of figure [3] We do not 
show results for data set C, as the number of SNe at z < 0.2 of that data set turned out to be insufficient for this type of test. 



In order to test the robustness of the fits we vary the value 
of the parameters Ay, cr v , cr,„, and the range in redshift (see 
table The most significant effect comes from a restriction 
of the redshift range to 0.02 < z < 0.2 or even 0.01 < z < 
0.1, which leads to a noticeable change in the best-fit value for 
qo and a substantial increase in its error. Shifting the maximal 
redshift to z — 0.3 changes the result only slightly. From this 
point of view it is reasonable to include all SNe with A v < 1 up 
to z — 0.2 in order to fit the Hubble law up to the second order 
in redshift. 

Although the full-sky fits indicate a robust data set C, our 
studies have shown that it is not well suited to be used within 
this work due to its small number of SNe and its bad sky cov- 
erage. We therefore focus on the results of data sets A, B and 
D in the following chapters. 

4.1.2. Galactic and equatorial hemispheres 

We start to study the anisotropy of the Hubble diagram by 
comparing those hemispheres which have a natural origin such 
as the galactic hemispheres and the equatorial hemispheres. 
The galaxy defines the zone of avoidance and is approxi- 
mately North-South symmetric. Thus the galactic hemispheres 
cover approximately the same observable regions of the sky 
and should be well suited for a cosmological test. On the 
other hand, most observations (especially at small redshifts) 
are done from the ground and thus the equatorial system is dis- 
tinguished. Any correlation with this system would hint to a 
systematic effect in the search, observation or data analysis of 
SNe la. 

Figure |4] shows the lcr contours for the galactic and equa- 
torial hemispheres. As expected, all three data sets show a con- 
sistent fit for North and South galactic hemispheres. However, 
the equatorial system shows some unexpected deviations. For 
data set A we find a significant deviation of the calibration 
Hu/Hq for the equatorial hemispheres. The northern hemi- 
sphere favours a lower value for Ho/ HI. The same trend is ob- 
served in data set B, but from the analysis of data set B alone we 



would not pay attention. For data set D the deviation is again 
clearly seen and is not restricted to the calibration. The corre- 
sponding best fit values and their errors are displayed in table[3] 
In order to quantify the evidence for an equatorial North- 
South asymmetry of the SN calibration, we use Monte Carlo 
(MC) simulations to check for artefacts of sky coverage. We 
test against 500 random realisations of the same data sets in 
which we mix the coordinates 1 and b for all SN la at z < 0.2. 
For each simulated hemisphere pair we calculate the deviation 
in^ 2 as 

A X 2 = mm(\x 2 N (S) - xl(N)l \x\{N) -^(S)|), 

where xjj(S ) denotes the value of expression (0 for the best- 
fit parameters of the North hemisphere applied on the data of 
the South hemisphere, etc. We also calculate the deviation in 
Ho I HI and qo of the best-fit values from the opposite hemi- 
spheres. The results of these hemispherical fits are given in ta- 
ble m 

While no deviation from isotropy is found for data set B, we 
find data set A shows a statistically significant anisotropy in the 
equatorial system. Only 0.8% of our MCs give rise to a larger 
difference in A% 2 and only 4.6% of the MCs show a larger dif- 
ference in the calibration Ho/Hq. The latter number for data set 
D is 5.0%. Thus, it seems the equatorial hemispheres do not 
agree with the expectation from the cosmological principle at 
the 95% C.L. with respect to the SN calibration (evidence from 
sets A and D). 

4.1 .3. Hemispheres of maximal asymmetry 

We found evidence for a significant asymmetry in the equato- 
rial coordinates in data sets A and D (and consistent with data 
set B). We may ask if the asymmetry is maximised by the equa- 
torial system or if more asymmetric directions exist on the sky. 
In order to test this anisotropy for its size we are going to iden- 
tify those hemispheres which have the largest deviation in^ 2 . 

We search for the maximal asymmetric hemispheres by 
rotating the poles across the sky. We calculate the deviation 
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in x 1 for each hemisphere with poles at I e [0°, 180°] and 
b e [-90°, 90°] in 1° steps. 

The results are given in tableland the deviations in Ay 2 are 
illustrated in Fig. [5] Each coloured field in the figure is 10°x 10° 
in size and its colour represents the mean deviation of the 100 
corresponding pole positions. It is striking that all three data 
sets give rise to the same pattern on the sky (especially sets 
A and D), despite the fact that the methods of SN data reduc- 
tion and selection criteria differ. By construction the pattern is 
symmetric and thus all information is contained in one hemi- 
sphere. In our discussion below we refer to the North galactic 
hemisphere. 

The asymmetry is strongest in data set A, followed by D 
and B. As data set A contains more SNe than set D, which con- 
tains more than set B, it seems that this trend is in accordance 
with the statistical power of the data sets. One can see, that two 
asymmetric directions are common in all three data sets. The 
first one is close to the equatorial poles at (/, b) = (123°, 27°), 
the other one is close to (/, b) = (10°, 15°). The proximity to 
the equatorial system is most pronounced in set A. Sets B and 
D maximise the asymmetry toward a direction ~ (70°, 15°). 

Let us also note that the observed pattern is unexpected. 
For randomly distributed SNe with a gaussian scatter in mag- 
nitude, we would expect that the resulting pattern would show 
less structure. 

We performed 500 MC simulations in order to check if the 
observed amount of asymmetry is to be expected for a data set 
with this kind of sky coverage. As above we mix the coordi- 
nates of the SNe of the data set under consideration and search 
for the maximal asymmetry in each of the simulated data sets. 
In order to keep the computational effort to a minimum, we 
now use a step size of 5°. The asymmetry of the MC data sets 
is compared to the one from the original set. The results are 
given in table |U The larger step size in the search for the max- 
imal asymmetric direction explains why the directions of max- 
imal asymmetry differ from the ones in table [3] For data sets 
A and D the asymmetry is larger than in 99.8% of our MCs. 
The deviation in the SN calibration seems to be significant at 
> 96% C.L. Data set D shows also a significant asymmetry in 
the extracted values of the deceleration parameter. Only 1.8% 
of the MCs show a larger difference. In contrast, the asymmetry 
in data set B appears not to be statistically significant. However, 
we should keep in mind that this data set shows a very similar 
asymmetry pattern but contains less SNe than sets A and D. 

4.1 .4. Summary of model-independent test 

For data set A we have identified a statistically significant 
asymmetry in the equatorial system. It turns out to be due to 
an off-set in the calibration Hq/H1 among the two hemispheres. 
The direction of maximal asymmetry is very close to that direc- 
tion in data set A. Data set B does not show the same amount 
of asymmetry, but is qualitatively consistent with data set A. 
However, it contains less SNe and it is thus expected that the 
asymmetry should be less obvious. The larger data set D shows 
again a statistically significant asymmetry in the SN calibra- 
tion. 




Fig. 5. Hemispherical asymmetry in A% 2 for three SN la data 
sets at small redshift (z < 0.2) for the model independent fit. 
When the pole lies in a region with a bluish colour the asymme- 
try is small, the red spots denote the directions of large asym- 
metry. We plot directions in galactic coordinates. The white 
points indicate the zenit (equatorial North) and its antipode 
(equatorial South). 



We cannot offer an explanation for the asymmetry, but we 
think that the proximity of the direction of maximal asymmetry 
to the equatorial poles suggests an systematic error in one of 
the steps (search, observation, data analysis) to the calibration 
of SNe. 

A second direction of asymmetry has been identified in all 
three data sets. Its detection is statistically less significant. 

4.2. Isotropy in the flat ACDM model 

Let us now include SNe la at any redshift in our test. This al- 
lows us to the test the concordance model of cosmology. We 
restrict our analysis to the flat ACDM model. Thus, the Hubble 
diagram can now be used to fit Qm and Hq/HL the calibration 
off-set of SNe la. Here we do not make use of data set D, as it 
is limited to small redshifts. 
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Fig. 6. Confidence contours for a full-sky fit to the flat ACDM 
model for three SNe la data sets. All SNe with Ay < 1 are 
included in the fits. We determine the calibration off-set Hq/Hq 
and the dimensionless mass density Qm- For comparison we 
show the corresponding contours of the WMAP measurement 
(Spergel et al. 2006). For the latter Hq indeed is the extracted 
Hubble parameter. Here the best-fit value of the calibration off- 
set is of no relevance, the vertical extension of the contours 
however reflects the physical scatter of the data. 
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4.2.1. Full sky 

To start with, we show the results for the full-sky fit in table 
[5] Again we see that the fits are quite robust, except for the 
restriction to z < 0.2. In that case, a matterless Universe (Q M = 
0) cannot be ruled out. The confidence contours for the ACDM 
fits to SNe la at arbitrary redshift are shown in figure [6] Only 
for data set C the best fit value is close to the one of the WMAP 
measurement. This is due to the fact that the SN calibration in 
set C was chosen to agree with WMAP measurements, whereas 
for SNe sets A and B no information from the CMB was used. 
In the following we restrict our presentation to data sets A and 
B, as the pencil beam geometry of data set C turned out not to 
be suitable for our tests. 



4.2.2. Galactic and equatorial hemispheres 

The fits to North and South galactic and equatorial hemispheres 
are presented in table [6] The asymmetry in equatorial coordi- 
nates found at low z is confirmed by the full data sets A and B. 
Again we run 500 MCs to test the statistical significance and 
find that the calibration off-set in the equatorial coordinates is 
now significant at > 95% C.L. in both data sets. In data set A 
it is at > 99% C.L. The corresponding confidence contours are 
shown in figure [7] 

Additionally, we find now that data set A shows also asym- 
metry in galactic coordinates, which is an asymmetry in the ex- 
tracted matter density at 95% C.L.; see table[6]for the strongly 
deviant values of in the North and South hemispheres, 
|AQmI = 0.2. This asymmetry is not at an significant level in 
data set B, but as for the model-independent test, data set B is 
fully consistent with the conclusions from data set A. 



Fig. 7. North (full lines) and South (dashed lines) confidence 
contours and best-fit values for galactic, equatorial and max- 
imum asymmetry hemispheres for the ACDM fit. These fits 
should be compared to the full-sky fits of figure [6] We do not 
show results for data set C, as the pencil beam geometry of that 
data set is not suitable for our test. 



4.2.3. Hemispheres of maximal asymmetry 

Finally, we search again for the maximally asymmetric pair of 
hemispheres. In figure|8]we show the hemispherical asymmetry 
Ax 2 in galactic coordinates. For data set A, the location of the 
maxima is very similar to the location of maxima in figure [5] 
and we confirm our findings above. 

In contrast, data set B shows a new structure, with a signifi- 
cant maximum close to (/, b) ~ (235°, 15°). This direction does 
not reflect any obvious large scale structure, but as we average 
over half of the sky, the pointing is not expected to be precise. 
However, an inspection of the asymmetry in the number of de- 
grees of freedom (dof) of the hemispheres shows that for data 
set B the maximum asymmetry in x 1 coincides with the maxi- 
mal asymmetry in the number of dof. We conclude that the pat- 
tern for data set B is due to the geometry of the SN sample and 
is fully consistent with an isotropic Hubble diagram. However, 
this finding does not rule out the precence of the effect observed 
in data set A for two reasons: data set A contains significantly 
more SNe and is dominated by nerby SNe (I = 0.27), whereas 
set B is dominated by objects at higher redshifts (I = 0.40). 

As above we run MCs to quantify the statistical significance 
of our findings. It turns out that the pattern and amount of asym- 
metry in set B is not unexpected, still data set A is unexpected 
at the 95% C.L. (see table©. 
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hemispherical 



data set A) 




Fig. 8. Hemispherical asymmetry in Ax 2 (see also figure[5]l in 
the context of the flat ACDM model. All SNe with A v < 1.0 at 
any redshift are included in the fits. 

4.2.4. Summary of the ACDM test 

We confirm the findings of the model-independent test, namely 
a significant asymmetry in the SN calibration (Ho/H^) between 
North and South equatorial hemispheres. On top of that asym- 
metries show up (especially in set B) that seem to be linked 
to the geometry of the samples. We do not find significant 
evidence for large scale structure effects at the largest scales 
(180°), but as objects like the Shapely supercluster are signifi- 
cantly smaller, our test is not well suited to identify such struc- 
tures. 

5. Results 

5.1. Is the Hubble diagram isotropic? 

We argue above that we find statistically significant deviations 
of Hubble diagrams form isotropy. These are most significant 
in data sets A and D, which contain most SNe. These asym- 
metries cannot be explained by the peculiar motion of the ob- 
server, as our analysis is done in the CMB rest-frame. If the 
asymmetry would be caused by peculiar motions of the SNe 
hosts, we would expect that its significance decreases when we 
include SNe at higher redshifts. This is not the case. As we go 
from the model-independent test at small redshifts (z < 0.2) 
to the model dependent test for arbitrary redshift, the signif- 
icance of the effect is actually increased. On top of that, the 
magnitude of the observed effect (10% in distance scale) ex- 
ceeds our expectations for the effect from peculiar velocities 
(v/(cz) ~ 10~ 3 /I), which in the most extreme case could ex- 
plain effects at the few per cent level in distance scale (or 
the SN calibration). However, it is not excluded that large 



scale bulk motions contribute to the effect, see e.g. Cooray & 
Caldwell (2006) who argue that a 5% variation of H might 
be possible in a low density bubble. Thus the origin of the 
anisotropy is most likely either due to a systematic effect in 
SN search, observation or data analysis or a (large) statistical 
fluke with a chance below 1:100. 

Evidence that we can also exclude the possibility of a sta- 
tistical fluke comes from the fact that the hemispheres of max- 
imal deviation have poles close to the poles of the equatorial 
system. This suggests that the origin of the observed anisotropy 
is, at least partly, due to systematic errors. A possible candidate 
for a systematic error would be inhomogeneous covering of the 
North and South sky, but our MC studies reveal that the asym- 
metry in the number of the objects cannot be held responsible. 
Exchanging the coordinates of the SNe (which preserves the 
asymmetry in number of objects) typically produces skies that 
are in agreement with the expectation of isotropy. Another pos- 
sible explanation might be that telescopes or search strategies 
in the North and South equatorial hemispheres have a system- 
atic calibration off-set. 

Yet another hint in favour of an unknown systematic ef- 
fect is the pattern observed in figures [5] and [8] If the Hubble 
diagrams would be isotropic, we would expect that there are 
several minima and maxima distributed randomly on the sky. 
However we observe quite regular patterns in the asymmetry 
maps. As the asymmetries in data set B are smaller than in data 
sets A and D, it seems better suited for cosmological analy- 
sis. The magnitude of its asymmetries could be consistent with 
large scale structure (1% effects), but are not statistically sig- 
nificant. However, the reason for the consistency of data set 
B with isotropy might just be that it does not contain enough 
SNe. Note that the fit values and their error bars of set B are 
consistent with the fits of the asymmetric sets A and D. 

Putting together all the evidence on the local flow direction 
in the Universe, mainly based on the analysis of galaxy clus- 
ters, Hudson et al. (2004) suggest a bulk flow of 225 km/s to- 
ward (/, b) = (300°, 10°) at depth greater than 60h l Mpc. One 
would expect that this should be reflected by SNe. However, 
our maximal asymmetric direction is significantly off-set from 
the large scale galaxy cluster bulk flow. A closer inspection of 
figure |5]for sets A and D shows a modest, statistically insignif- 
icant, asymmetry close to the suggested bulk flow direction. 

5.2. Is the local Hubble expansion rate different from 
the global one? 

We cannot answer that question, as we have no handle in our 
analysis to the absolute calibration of SNe. If one could exclude 
that the calibration asymmetries found in this work are due to 
systematic errors, this would be a strong indication for a vari- 
ation of the Hubble rate of the order of 10%. If so, the tension 
between the WMAP 3yr measurment of Hq and the determina- 
tion based on SNe by Sandage et al. (2006) might be reconciled 
with each other. 

The order of magnitude would also be consistent with 
recent claims by Jha et al. d20061 l about the existence of 
a "Hubble bubble"; but more recent analysis (Wang 120071 
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Conley et al. 2007) does not confirm the claim. Conley et 
al. (|2007) show that the evidence for a "Hubble bubble" de- 
pends crucially on how SN colours are modelled. As we think 
that our test points towards an equatioral North/South system- 
atic, we do not think that this work should be regarded as a 
support of the idea of a local bubble. However, it might be 
interesting to check if an equatorial systematic in SN colours 
could be found. 

5.3. Is the Universe accelerating today? 

Based on our analysis one cannot answer this question with a 
straight 'yes' . The model-independent analysis is fully consis- 
tent with an accelerated Universe, but the evidence in favour of 
acceleration is at most at the 2<x level. For the supposedly most 
accurate compilation of local SNe from Jha et al. 2006 (data set 
D), we find that qo — is within the 2cr contour (see our figure 
0. The same is true for all four analysed sets of data (includ- 
ing the SNLS data set, not shown in our figure [3). What can 
be confirmed is that the Einstein-de Sitter model (go = 1/2) is 
outside the 2cr region, and is therefore disfavoured already by 
the model-independent test at small redshift. 

If we assume that the flat ACDM is the full truth, then the 
SNe la data can indeed provide convincing evidence that the 
Universe is accelerating. Within that model deceleration (qo > 
0) corresponds to Q.m > 2/3. As can be seen from figure [6] this 
case is excluded at high confidence. However, we should keep 
in mind that neither the physics of 'A', nor that of 'CDM' is 
understood at any depth. 

6. Conclusion 

The purpose of the presented study is to develop tools for 
and to test the isotropy of Hubble diagrams, which are at the 
very foundation of modern cosmology. Within our established 
cosmological model we expect some small deviations from 
the isotropy, but with the present day accuracy of SNe ob- 
servations, we would not expect to be able to detect them at 
high statistical significance (apart from our proper motion). 
Nevertheless, we set out to apply the hemisphere tests to ex- 
isting data in order to develop the methodology and to test for 
systematic effects. 

To our surprise we identified a statistically significant 
asymmetry, which is maximised close to the orientation of the 
equatorial system. It seems to us that this calls for a thorough 
investigation of possible systematic effects, which is beyond 
the scope of this work. Our analysis indicates that there is an 
off-set in the calibration of SNe between the equatorial hemi- 
spheres. Typically SN searches are flux-limited (not red-shift 
limited) and thus one could imagine that the search in the North 
and in the South selects samples with different dispersion and 
mean value in SN brightness. Our findings would be consistent 
with a more complete search for nearby SN in the South com- 
pared to the North. We think that this study shows, besides the 
interest in the large scale structure, that a large sky coverage of 
SN search missions is also important for the issue of systematic 
errors. 



The Hubble diagram is currently the only direct mean to 
probe the acceleration of the Universe. Most of our evidence 
for the present day acceleration comes from indirect arguments 
and relies on a bunch of untested assumptions. Our model- 
independent test fails to detect acceleration of the Universe at 
high statistical significance. It seems to us that it is too early 
to take accelerated expansion of the Universe for granted, as 
the evidence heavily relies on the a priori assumption of the 
ACDM model. 

Acknowledgements. We thank Camille Bonvin, Ruth Durrer, Chris 
Gordon, Steen Hannestad, Stefan Hofmann, Michael Hudson, Ariel 
Goobar, Martin Kunz, Bruno Leibundgut and Anze Slozar for discus- 
sions and comments. 

References 

Astier, P., et al. 2006, A&A, 447, 31 larXiv:astro-ph/05 10447 1 
Barris, B.J., et al. 2004, ApJ, 602, 571 |arXiv:astro-ph/03 10843 | 
Bennett, C.L., et al. 2003, ApJS, 148, 1 |arXiv:astro-ph/0302207 1 
Bonvin, C, Durrer, R., & Gasparini, M.A. 2006a, Phys. Rev. D, 73, 

023523 |arXiv:astro-ph/0511183| 
Bonvin, C, Durrer, R., & Kunz, M. 2006b, Phys. Rev. Lett., 96, 

191302 |arXiv:astro-ph/0603240| 
Conley, A., et al. 2007, ApJ, 664, L13 larXiv:0705.0367l 
Cooray, A., & Caldwell, R.R. 2006, Phys. Rev. D, 73, 103002 

|arXiv:astro-ph/0601377) 
Copi, C.J., Huterer, D., Schwarz, D.J., & Starkman, GD. 2007, Phys. 

Rev. D, 75, 023507 |arXiv:astro-ph/0605 1 35 1 . 
Eriksen, H.K., Hansen, F.K., Banday, A.J., Gorski, K.M., & 

Lilje, P.B. 2004, ApJ, 605, 14 [Erratum-ibid. 609, 1198] 

|arXiv:astro-ph/0307507l . 
Freedman, W.L., et al. 2001, ApJ, 553, 47 |arXiv:astro-ph/0012376| . 
Germany, L.M. 2001, Ph.D. Thesis, Australian National University 
Gordon, C, Land, K., & Slosar, A. 2007, larXiv:0705.1718l 
Gupta, S., Saini, T.D., & Laskar, T. 2007, |arXiv:astro-ph/0701683l 
Guy, J., Astier, P., Nobili, S., Regnault, N., & Pain, R. 2005, 

|arXiv:astro^ph /Q506583 
Hannestad, S., Haugb0lle, T., & Thomsen, B., 2007, larXiv:0705.0979l 

[astro-ph]. 

Haugb0lle, T, et al. 2006, |arXiv:astro-ph/0612"T37l 

Hudson, M.J., Smith, R.J., Lucey, J.R., & Branchini, E. 2004, 

MNRAS, 352, 61 |arXiv:astro-ph/0404386) 
Hui, L., & Greene, P.B. 2006 Phys. Rev. D, 73, 123526 

|arXiv:astro-ph/0512l59"l . 
Jha, S. 2002, Ph.D. Thesis, Harvard University 
Jha, S., Riess, A.G., & Kirshner, R.P 2006, ApJ, 659, 122 

|arXiv:astro^ph /0612666 1 
Kolatt, T.S., & Lahav, O. 2001, MNRAS,323, 859 

|arXiv:astro-ph/000804T) 
Neill, J.D., Hudson, M.J. & Conley, A., 2007, larXiv:0704.1654i rastro- 

ph]. 

Perlmutter, S„ et al. 1999, ApJ, 517, 565 | arXiv:astro-ph/9812133 | 
Phillips, M.M., et al. 1999, AJ, 118, 1766 |arXiv:astro-ph/9907052 1 
Radburn-Smith, D.J., Lucey, J.R. & Hudson, M.J., 2004, MNRAS, 

355, 1378 |arXiv:astro-ph/040955 1 1 
Riess, A.G., Press, W.H. & Kirshner, R.P. 1995, ApJ, 445, L91 

|arXiv:astro-ph/9412017l 
Riess, A.G., Davis, M., Baker, J. & Kirshner, R.P. 1997, ApJ, 488, LI 

|arXiv:astro-ph/970726T) 
Riess, A.G., et al. 1998, AJ, 116, 1009 farXiv:astro-ph/9805201 J 
Riess, A.G, et al. 2004, ApJ, 607, 665 |arXiv:astro-ph/0402512 | 



Dominik J. Schwarz and Bastian Weinhorst: (An)isotropy of the Hubble diagram: comparing hemispheres 1 1 





dof 


2 
r 

dof 


K 


9o 


data <ipt A - 753 SNe 7 F TO 00? 1 7551 










4.. < 1 rr - km/s 7 < ? 


137 


1.27 


1 0? + 0? 


-0 78 + 90 


Ay < 0.5 


116 


1.30 


1 0? + 0? 


-0 56 + 94 


rr - 730 km/s 

iy y — z 'V; iviii/ a 


137 


1.82 


1 03 + 0? 

i.uj in v/.v/z. 


_0 Q7 + 85 


rr — 4fS0 Ifm/s 

ly y — T^VJVy Ivlll/ a 


137 


0.97 


1 0? + 03 


-0 68 + 96 


(J y — U7U Mil/ a 


137 


0.67 


1 0? + 03 


-0 61+1 06 


07 < 7 < 7 


73 


1.20 


1 01 + 03 


-0 3 1 + 1 08 


z < 0.1 


128 


1.33 


1.02 ± 0.03 


-0.52 ± 0.52 


01 < 7 < 1 


97 


1.11 


1 00 + 03 


0.28 ± 1.25 


z < 0.3 


142 


1.23 


1.02 ± 0.02 


-0.65 ± 0.73 


data <ipt R- 1 86 SNe 7 P TO 01 1 7551 










A < 1 /-T- _ -245 km/s 7 < ? 

./t. \/ _ J. ; L» y — jtj ivlll/ a, yj.z. 


75 


0.84 


1 01 + 03 


-1.42 + 1.23 


A v < 0.5 


66 


0.79 


1 00 + 04 


-1.36 + 1.79 


rr - 730 km/s 


75 


0.92 


1 01 + 03 

1.U1 ^ v/. v/ J 


-1.43 + 1.20 


rr — 4(S0 km/s 

U y — T^VJVy Ivlll/ a 


75 


0.75 


1 01 + 03 

L,\J± Vy . Vy J 


-1.40 + 1.25 


n- — AQO km/s 

(7 i' — U7U Mil/ a 


75 


0.60 


1 01 + 04 

1 . v/ 1 ZC Vy . Vy^ 


_1 -ic . 1 n 

1 . *J ZLL 1 . J»y 


0.02 < z < 0.2 


50 


0.95 


1.01 ± 0.04 


-1.27 ± 1.36 


z < 0.1 


70 


0.84 


1 02 + 05 


-2.16 ± 2.21 


z < 0.3 


80 


0.81 


1 00 + 03 

i » vy vy _!_ vy . vy ^y 


-0 78 + 89 


data «spt C- 1 17 SNe rf TO 01 5 1 011 

U 11 Kl 3CL V-- * 11 / J1>C, 4j vZ [U.U1 J, 1 . V7 1 J 










(7 y — RJ.11/a, u mi — u.UJ) 4, V/.z. 


42 


0.84 


1 07 + 04 

/ m Vy . Vy ^ 


-0 57 + 1 63 

Vy . J / I 1 .UJ 


rr - 230 km/s 

17 y — Z-^y Vy Ivlll/ a 


42 


0.92 


1 08 + 04 

1 .\J(J _!_ Vy . Vy^ 


-0.57 ± 1.59 


rr — 4rtD km/s 
17 y — t^uu iviii/ a 


42 


0.76 


1 08 + 05 

i.uo ^ vj.vj»y 


-0.58 + 1.71 


rr — 600 km/s 

(J y — U7U Kill/ a 


42 


(SO 


1 08 + D Dri 
1 .V70 i u.uu 


-0 (S^ + 1 88 


0? < 7 < 7 


30 


1.07 


1 .\J 1 ~J1 Vy . Vy^/ 


-0 78 + 1 83 

\J. Z.O i 1 . J 


7 < 1 


40 


0.87 


1 06 + 05 


16 + 7 15 


z < 0.3 


46 


0.80 


1.08 ± 0.03 


-0.56 ±0.58 


rr , = 

u ml u 


42 


3.73 


1 09 + 03 

1 ivy —i— vy > vy ^y 


-0.82 ± 0.75 


rr- . - 02 


42 


1.40 


1 08 + 04 

J. .v/(j -!- Vy . Vy^ 


-0 62 + 1 27 

Vy . v." i_ _L_ 1 . / 


data «spt 1 31 SNe 7 <= TO 00? 1741 










A„ < 1 rr - 345 km/s rr -0 016 7 < 7 

Jl y _Z J. • L* \' — JT^J Mil/ a, 17 — Vy . Vy 1 17, _ \/._- 


117 


1.37 


1 01 + 03 

l.V/1 Vy . Vy ^7 


_i 39 + 1 35 


A v < 0.5 


99 


1.43 


1 00 + 0? 

1 . VJvy ~J1 Vj.VjZ. 


-0 86 + 1 36 


rr - 730 km/s 

17 y — Z.JU iviii/ a 


117 


2.14 


1 0? + 0? 
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rr — 4rtf) km/s 

17 y — T^VJVy IV111/ a 


1 17 


0.98 
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1.U1 ^ v/. v/ J 


-1.16 + 1.45 
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0.59 
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-0 86 + 1 60 


0? < 7 < 7 

U.UZi V, ,, 17. Z, 


50 


1.24 


Q7 + 04 


0.11 + 1.59 


7 < 1 
<, _i V7. I 


1 15 


1.39 


1 0? + 03 


-1.49 + 1.71 


0.01 < z < 0.1 


84 


1.08 


1.00 ±0.03 


-0.81 ± 1.78 


Cr,,,/ = 


117 


1.64 


1.01 ±0.02 


-1.30 ± 1.20 


cr,,,, = 0.03 


117 


1.05 


1.02 ±0.03 


-1.51 ± 1.66 



Table 2. Robustness of the full-sky fit of the calibration Hq / H* () and the deceleration parameter qo at small redshifts. We compare 
the number of degrees of freedom (dof), ^- 2 /dof, and the best fit cosmological parameters for the four data sets described in the 
text for various assumptions on acceptable light extinction Av, peculiar velocity dispersion cr v , intrinsic dispersion cr mt , as well 
as redshift interval included in the fit. Our analysis with Ay < 1 includes all SNe without information on Ay, but we exclude 
those when investigating Ay < 0.5. 
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Table 3. Hemisphere fit of Hq/H^ and qo at small redshifts. We compare the number of degrees of freedom (dof), v; 2 /dof, and the 
best fit cosmological parameters for the data sets A, B and D for the galactic, equatorial and maximal asymmetry hemispheres. 
In addition we show for the galactic and equatorial hemispheres the percentage of MC simulations with a larger deviation in^- 2 , 
Hq/Hq and qo. In brackets we provide the respective differences Ay 2 , and Ago- 
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(0.13) 


(5.62) 



Table 4. Hemispheres of maximal asymmetry from a 5° search (in contrast to the more accurate 1° search in tableO. The purpose 
of the sparse search is to run MC simulations. The meaning of the columns is explained in the caption of table[3] 
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data set A: 253 SNe, z e [0.002, 1.755] 
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data set B: 186 SNe, z e [0.010, 1.755] 
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42 


0.84 


1.08 ±0.04 


20+ 148 


f&a = 


115 
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Table 5. Robustness of the full-sky fit to the flat ACDM model at arbitrary redshift. The fit parameters are the calibration 
Hq/Hq and the dimensionless matter density £2m- We compare the number of degrees of freedom (dof), ^ 2 /dof, and the best fit 
cosmological parameters for data sets A, B, and C for various assumptions on the acceptable light extinction Ay, peculiar velocity 
dispersion cr v , intrinsic dispersion cr int and redshift interval included in the fit. 
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Table 6. Comparison of galactic, equatorial and maximal asymmetry hemispheres for the flat ACDM model. We use all SNe 
with Ay < 1.0 from data sets A and B. The meaning of the columns is explained in the caption of table [3] 



